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ABSTRACT 

Using RHESSI hard X-ray imaging spectroscopy observations, we analyze 
electron flux maps for a number of extended coronal loop flares. For each event, 
we fit a collisional model with an extended acceleration region to the observed 
variation of loop length with electron energy E, resulting in estimates of the 
plasma density in, and longitudinal extent of, the acceleration region. These 
quantities in turn allow inference of the number of particles within the acceler- 
ation region and hence the filling factor / - the ratio of the emitting volume to 
the volume that encompasses the emitting region(s). We obtain values of / that 
lie mostly between 0.1 and 1.0; the (geometric) mean value is / = 0.20 x/-r 3.9, 
somewhat less than, but nevertheless consistent with, unity. Further, coupling 
information on the number of particles in the acceleration region with informa- 
tion on the total rate of acceleration of particles above a certain reference energy 
(obtained from spatially-integrated hard X-ray data) also allows inference of the 
specific acceleration rate (electron s _1 per ambient electron above the chosen ref- 
erence energy). We obtain a (geometric) mean value of the specific acceleration 
rate 77(20 keV) = (6.0 x /-=- 3.4) x 10~ 3 electrons s" 1 per ambient electron; this 
value has implications both for the global electrodynamics associated with replen- 
ishment of the acceleration region and for the nature of the particle acceleration 
process. 
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gamma-rays 



1 Dipartimcnto di Matematica, Universita di Genova, via Dodecanese) 35, 16146 Genova, Italy; 
guo@dima. unige . it , piana@dima. unige . it 

2 Department of Physics and Astronomy, Western Kentucky University, Bowling Green, KY 42101; em- 
slieg@wku.edu 

3 CNR - SPIN, via Dodecaneso 33, 1-16146 Genova, Italy; annamaria.massone@cnr.it 



-2 - 



Introduction 



An important diagnostic of high-energy electrons accelerated in solar flares is the hard 
X-ray bremsstrahlung that they produce as they propagate through the ambient solar atmo- 
sphere. The Ramaty High Energy Solar Spectroscopic Imager (RHESSI ) has revealed a new 
class of flares in which the bulk of the hard X-ray emission is prod uced predominantly not 
in dense chromospheric footprints, but rather in the coronal loop ( jVeronig fc Brown! 12004 ; 
Sui et al.l |2004i iKrucker et al.l 120081 ). For such sources, the corona is not only the site of 
particle acceleration, but also dense enough to act as a thick target, stopping the accelerated 
electrons before they can penetrate to the chromosphere. 

For suprathermal electrons with energy substantially greater than the thermal energy 
of the ambient electrons with which they in teract, it is appropriate to use a collisional cold- 
target energy loss rat e (e.g., Emsliel 1978 ). for which the penetration depth of electrons 



increases with energy. IXu et al.l (l2008h analyzed a set of extended coronal flare loops located 
near the solar limb, and were indeed able to account for the observed behavior of loop extent 
with ph oton energy e in t erms of a cold-target collisional model with an extended acceleration 



region. iGuo et al.l ( 120121 ) have extended this analysis technique to a study of the variation 
of loop size with electron energy E, in which the visibilities used to construct the electron 
flux ima ges are obtained b y regularized spectral inversion of the visibility data in the count 
domain (jPiana et al.l 120071 ) . 



Here we apply this new analysis technique to several simple coronal loop events ob- 
served by RHESSI. In Section [21 we present basic data for the 22 events used in the study. 
In Section B] we fit t he variation of loop size with electron energy E to the parametric model 
of IGuo et al.l (120 121 ) in order to determine the acceleration region length L and density n 
for each event. In Section H] these values are used to determine estimates of two important 
properties of the acceleration region - the filling factor / (the ratio of the volume that is 
actively involved in electron acceleration to the overall volume that encompasses the acceler- 
ation region [s]) and the specific acceleration rate t](Eq) (the rate of acceleration of electrons 
to energies > E per ambient electron), and we compare the values of these quantities to the 
predictions of various acceleration models. 



2. Events Studied 



The list of events studied is shown in Table [TJ In this context, an "event" is a time 
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Table 1. Event List and Spectral Fit Parameters 



Event No. 


Date 


Time (UT) 


EM (10 49 cm" 3 ) 


T (keV) 


S 


E t (keV) 


dM/dt (10 35 s- 1 ) 


1 


2002-04-12 


17:42:00-17:44:32 


0.30 


1.53 


8.24 


15.5 


2.71 


2 




17:45:32-17:48:00 


0.46 


1.54 


8.01 


15.5 


4.69 


3 


2002-04-15 


00:00:00-00:05:00 


0.22 


1.75 


7.48 


15.5 


4.70 


4 




00:05:00-00:10:00 


0.76 


1.61 


7.93 


15.5 


9.32 


5 




00:10:00-00:15:00 


1.02 


1.60 


8.37 


15.5 


11.41 


6 


2002-04-17 


16:54:00-16:56:00 


0.06 


1.51 


5.70 


15.5 


0.39 


7 




16:56:00-16:58:00 


0.22 


1.43 


8.78 


14.8 


2.43 


8 


2003-06-17 


22:46:00-22:48:00 


1.92 


1.71 


9.95 


16.5 


17.27 


9 




22:48:00-22:50:00 


2.59 


1.67 


10.36 


16.5 


17.91 


10 


2003-07-10 


14:14:00-14:16:00 


1.26 


1.45 


10.05 


15.5 


7.43 


11 




14:16:00-14:18:00 


1.31 


1.34 


10.38 


14.8 


8.53 


12 


2004-05-21 


23:47:00-23:50:00 


0.35 


1.85 


7.07 


18.5 


3.28 


13 




23:50:00-23:53:00 


0.62 


1.75 


7.51 


18.5 


2.32 


14 


2004-08-31 


05:31:00-05:33:00 


0.06 


1.61 


10.56 


15.5 


0.40 


15 




05:33:00-05:35:00 


0.21 


1.57 


12.19 


18.5 


0.29 


16 




05:35:00-05:37:00 


0.29 


1.48 


7.45 


18.5 


0.16 


17 


2005-06-01 


02:40:20-02:42:00 


0.14 


1.81 


6.53 


17.5 


1.44 


18 




02:42:00-02:44:00 


0.37 


1.70 


7.86 


17.5 


2.67 


19 


2011-02-13 


17:33:00-17:34:00 


0.54 


1.39 


5.86 


10.5 


35.02 


20 




17:34:00-17:35:00 


0.52 


1.68 


6.55 


14.5 


19.43 


21 


2011-08-03 


04:31:12-04:33:00 


0.36 


1.61 


9.23 


15.5 


3.96 


22 


2011-09-25 


03:30:36-03:32:00 


0.13 


1.44 


8.33 


14.5 


1.19 
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Fig. 1. — Left panel: Light curves, in the energy intervals labeled at the top right of the plot, 
for the flare on 2002 April 17. The vertical lines delineate the time intervals for Events 6 
and 7. Right panel: Spectral fit to Event #7 (16:56:00 - 16:58:00 UT). The green histogram 
shows the thermal component of the spectrum (EM = 0.216 x 10 49 cm -3 ; T = 1.43 keV) 
and the yellow histogram shows the non-thermal thick-target component (transition energy 
E t = 14.8 keV; spectral index 5 = 8.78). The red histogram represents the sum of the 
thermal and nonthermal components and the lilac histogram represents the background. 
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interval during a flare for which spatial and spectral observations are sufficiently good to 
permit both a determination of the source spatial structure at a variety of energies and the 
overall spectrum of the hard X-ray emission. Some flares provide multiple "events"; other 
flares only one (see Table [I]). For each event, we fit the spatially-integrated hard X-ray 
emission with an isothermal-plus-power-law form, yielding values (Table [1]) of the emission 
measure EM (cm~ 3 ) and temperature T (keV) of the thermal source, the intensity and 
spectral index 6 = 7 + 1 of the injected nonthermal electron spectrum (corresponding to the 
hard X-ray spectral index 7), and E t (keV), the transition energ y between th e thermal and 



nonthermal components. Straightforward thick-target modeling (lBrownlll97l[ ) then provides 
dM/dt (s _1 ), the acceleration rate of electrons above (somewhat arbitrary) reference energy 
E = 20 keV. 

Parenthetically, we note that all the injected electron spectra are rather steep (the 
lowest value of 5 is 5.70 [Event #6], and the typical value is in the range 7-9). While 
this may indicate a property of the electron acceleration process in a relatively dense (see 
below) medium, it may also simply be an observational selection effect - the relative paucity 
of high-energy electrons in such steep spectra is consistent with the absence of footpoint 
emission that such high-energy electrons would produce. 

Figured] shows the RHESSIcowat rate profiles for Events 6 and 7 (16:54:00 - 16:56:00 UT, 
and 16:56:00 - 16:58:00 UT, respectively, on 2002 April 17) in five different energy channels. 
We have identified with vertical lines the time intervals for each event. The right panel shows 
the spectrum for Event #7, with the values of the spectral fit parameters provided in the 
caption (and in Table [p. 



Electron flux images of each event, in the representative 18-20 keV energy channel, are 
shown in Figure |2j These images were produced by performing a spectral inversion on the 
count visibility dat a to obtain the corre sponding electron visibilities and then using the uv- 
smooth algorithm ( Massone et al. 2009h to produ ce maps of the electron flux (weighted by 



the line-of-sight column density; iPiana et al. 



20071 ). Two aspects of this technique are worthy 



of note. First, the technique exploits the fact that for the bremsstrahlung process counts of 
energy q are produced by all electrons with energy E > q, so that the method yields images 
at electron energies E up to and beyond the maximum count energy g max observed. Second, 
the regularized spectral inversion procedure that produces the electron visibilities results in 
images that, by construction, vary more smoothly with electron energy E than the "parent" 
count images vary with q. This smooth variation with E greatly facilitates the analysis of 
the following section. 
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3. Analysis 



Using the electron flux maps, we first calculate the princi pal longitudinal a nd lateral 
directions s and t for the source, using the technique described in lGuo et all (120121 ) . The lon- 
gitudinal extent of the source external to the acceleration region can be found by considering 
the standard deviation 



a{E) 



if™s 2 F(E,s) ds 



f °°F(E,s)ds ' 

where F(E, s) is the electron flux spectrum at longitudinal position s. 

Various physical processes can in principle contrib ute to the be havior of F(E,s). Ob- 
viously, Coulomb collisions wit h ambient electrons (e.g., Emsliel 19781 ) must be considered. In 



addit ion, various authors (e.g.. lKnight fc Sturrocklll977 



Emsliel 



1981 



Zharkova fc Gordovskyy 



20061 ) have stressed the need to consider also the Ohmic energy losses associated with driving 
the beam-neutralizing return current in a resistive medium. Such energy losses are propor- 
tional to the decelerating voltage difference, which is in turn proportional to the beam 
current. Ret urn current losses are therefore likely to be significant only in large events (e.g., 
Emsliel 1 198 11 ) with electron acceleration rates dM/dt substantially greater than those con- 
sider ed here (Table ED) . For similar reasons, we be lieve that collective plasma effects (see, 



Hoyng fc Melrose 1977 : Emslie fc Smith 19841 ) are also likely to be unimportant. We 



e.g. 

therefore consider only Coulomb energy losses in the computation of F(E, s). 

We therefore consider a cold-target collisional injection model and a target of uniform 
density n (cm -3 ). We also neglect pitch angle scattering an d dispersion around mean values, 
which typically affect F(E, s) by factors of order unity (IBrownl Il972l ; iLeach fc Petrosian 
19811 ). and we shall address these factors briefly in Section HI For such a scenario, the form 
of F(E, s) can be deduced from the one- dimensional continuity and energy loss equations 



F(E) dE = F (E ) dE ; ^ = 

ds E 



(2) 



(Here K = 27re 4 A, e being the electronic charge and A being the Coulomb logarithm.) The 
solution for E(s) is E 2 = Eq — 2Kns, so that dE^/dE = E/Eq and hence, for a power-law 
injection spectrum F (E ) ~ Eq 5 , 



F(E,s) 



E 



(E 2 + 2Kns)( s + 1 )/ 2 ' 
Substituting the expression ([3]) into Equation (TTT) , we obtain, after some algebra, 



(3) 
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For a mode0in which electrons are accelerated within a region extending from [— L /2,L /2] 
and injected into an external region with uniform density n, we therefore arrive at the 
relationship between the observed longitudinal source extent L and electron energy E: 



L(E) 




E\ (5) 



2 2 KnV (<J-3)(<5-5) 

For each event, the values of the acceleration region length L and the loop density n are ob 
tained by best-fitting Equation (J5]) to the inferred form of the variation of loop le ngth L(E) 



with electron energy E in the predominantly nonthermal domain E>E t (see iGuo et al. 



20121 ) . The resulting best-fit parameters are presented in Table [2j From the electron flux 
images, we also straightforwardly obtain the width (lateral extent) W of the emitting re- 
gion, which typica lly exhibits a much smaller variation with energy E than does L (see 



Kontar et al.ll201ll ) and can hence be taken as a constant. From the inferred values of Lq, 



W, and n, we obtain the volume of the acceleration region 



Vo = — ^ — (6) 

and the number of particles it contains 

M = nV Q . (7) 
Values of Vq and M are provided for each event in Table [21 



3.1. Specific Acceleration Rate 



The specific acceleration rate (electrons s _1 per electron) is defined (jEmslie et al.ll2008l ) 
as the ratio of two quantities: dj\f/dt(> E Q ), the rate of acceleration of electrons beyond 
energy Eq, and A/", the number of particles available for acceleration: 



l - Xu et al. ( 20081) also discuss a more physically self-consistent model which incorporates the finite density 



in the acceleration region in the form for L(E). The corresponding expression for L(E) cannot be expressed as 
a simple closed form and hence is more complicated to use in a best-fit analysis. However, this more correct 
form nevertheless yields results for Lo and n that are comparable to those obtained from the "tenuous 
acceleration region" model used here. 



- 9 - 



The quantity dj\f/dt(> E ) is readily determined by spectral fitting of the spatially-integrated 
hard X-ray emission - see Table [H The quantity Af can be found from Equation ([7]) - see 
values in Table |2J We can thus deduce the value of the specific acceleration rate r](E ) in 
each event; values are given in Table [2 



3.2. Filling Factor 

The soft X-ray emission measure EM is related to the plasma density n and the emitting 
volume Kmit through EM = n 2 V emit . Given that an emitting region may be composed of 
a number of discrete emitting subregion (e.g., "strands," "kernels"), the emitting volume 
may be equal to or smaller than the total flare volume V estimated from observations of 
the spatial extent of the source. The ratio of the emitting volume to the volume V of the 
observed region that encompasses the emitting region(s) is termed the filling factor 

, _ Kmit _ EM 

V n V 

For each event, we estimated / by using the value of EM from the spectral fit to the thermal 
portion of the hard X-ray spectrum (Table [1]), the density n from the fit to Equation ([5]) 
(Table [2]), and the observed volume V of the emitting portion of the loop, measured at the 
transition energy E t (Table [TJ, the maximum energy at which thermal emission is predomi- 
nant: V = (vr/4) W 2 L(E t ). Values of / for each event are given in Table [2j 



4. Results and Conclusions 

The values of L , W, V , n, Af, 7/(20 keV) and / for each event are presented in Table [2J 
While statistical uncertainties in these values could readily be calculated through a Monte 
Carlo metho d in which noise is added to the RHESSI count visibility data and the process 



repeated (see lGuo et al.ll2012l ). we have intentionally refrained from doing so here, since the 
approximations and assumptions used in the model doubtless entail even larger uncertainties. 
Instead, we let the scatter of the inferred values of the parameters across the 22 events 
determine the extent over which the parameters range. We have calculated (Table [2]) the 
value of the (geometric) mean value of each quantity and the (multiplicative) uncertainty in 
this value. In particular, we obtain n = (2.9 x/-r 1.9) x 10 11 cm -3 , / = 0.20 x/-=- 3.9, and 
?7(20 keV) = (6.0 x /-r 3.4) x 10~ 3 electrons s _1 per ambient electron. 
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Table 2. Acceleration Region Characteristics 



Event No. 


Lo (arcsec) 


W (arcsec) 


Vo (100 arcsec 3 ) 


n (10 11 cm" 3 ) 


M (10 37 ) 


r, (20 keV) (10~ 3 s" 1 ) 


/ 


1 


18.6 


7.0 


7.2 


1.5 


4.1 


6.5 


0.45 


2 


16.3 


6.9 


6.2 


1.4 


3.2 


14.5 


0.83 


3 


16.7 


7.3 


7.0 


4.4 


11.7 


4.0 


0.04 


1 


16.6 


7.3 


7.0 


4.8 


12.8 


7.3 


0.11 


5 


16.6 


8.2 


8.7 


10.5 


34.9 


3.3 


0.03 


6 


11.9 


5.9 


3.3 


4.9 


6.0 


0.6 


0.02 


7 


10.4 


6.0 


3.0 


1.8 


2.0 


12.1 


0.44 


8 


17.8 


6.9 


6.4 


2.6 


7.1 


24.1 


0.90 


9 


18.8 


6.6 


6.5 


2.9 


7.7 


23.1 


1.05 


10 


15.1 


6.0 


4.2 


2.9 


5.4 


13.8 


0.72 


11 


16.0 


5.7 


4.1 


1.9 


3.1 


27.8 


1.95 


12 


10.3 


6.6 


3.5 


5.1 


6.7 


4.9 


0.08 


13 


9.9 


6.5 


3.3 


4.6 


5.7 


4.1 


0.18 


14 


21.5 


5.3 


4.8 


1.5 


2.8 


1.4 


0.13 


15 


17.4 


6.3 


5.4 


0.8 


1.7 


1.7 


1.03 


16 


17.8 


6.4 


5.8 


2.3 


5.1 


0.3 


0.18 


17 


11.0 


6.2 


3.3 


3.9 


5.0 


2.9 


0.05 


18 


9.9 


6.3 


3.1 


3.2 


3.8 


7.0 


0.22 


19 


19.9 


6.2 


6.1 


11.1 


25.7 


13.6 


0.02 


20 


14.5 


6.1 


4.2 


5.2 


8.3 


23.4 


0.10 


21 


9.9 


6.1 


2.9 


2.2 


2.4 


16.5 


0.53 


22 


12.4 


6.0 


3.6 


1.7 


2.3 


5.2 


0.26 


Geometric Mean 


14.5 


6.4 


4.7 


2.9 


5.4 


6.0 


0.20 


X/- 


1.3 


1.1 


1.4 


1.9 


2.2 


3.1 


3.9 
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Returning to the simplifying assumptions used in determining the form of the electron 
flux F(E, s) (Equations [2] and [3]), we note that inclusion of electron trajectories that have 
a non-zero pitch angle to the guiding magnetic field and/or a guiding magnetic field that is 
inclined to the longitudinal axis (the direction defining the coordinate s) will add a factor 
fi = cos 8, where 9 is the angle between the electron velocity vector and the longitudinal 
direction, to the energy-dependent term in Equation ([5]). This will result in a decrease (by a 
multiplicative factor /i) in the inferred density n, which in turn, by Equations ([7]), ([8]) and fl9]), 
will increase the values of / and 77 by factors of 1//J 2 and 1//j, respectively. Inclusion of return 
current Ohmic energy losses and/or energy losses to waves through collective plasma effects 
will also decrease the electron penetration depth, leading to further decreases in the inferred 
value of n and so increases in / and rj. The values of / and 77 cited above are therefore in 
all likelihood lower limits. 

The inferred values of / are generally somewhat less than unity, with the exception of 
three events (## 9, 11 and 15), for which / = 1.05, 1.95 and 1.03, respectively. Given the 
uncertainties in the data, the approximations in the analysis method, and the factor of four 
spread in the inferred values of /, neither of these values exceeds unity by an alarming margin. 
The mean value of the filling factor obtained is consistent, within a logarithmic standard 
deviation or so, with unity. This result, while not entirely surpris ing, is nevertheles s still 



significant. It validates the assumption used by many authors (e.g.. lEmslie et al.ll2004l ) that 
most of the observed flare volume contains bremsstrahlung-emitting electrons; the degree 
to which the emission is fragmented (e.g., striated into "kernels" or "strands" of emission 
situated within a relatively inert background medium) is quite small. 

The inferred mean value of 7/(20 keV) ~ 5 x 10 -3 electrons s _1 per ambient electron is 
broadly cons i stent with the values reported for a series of extended-loop-source events by 
Emslie et al.l (120081 ) . It should also be noted that the value of the specific acceleration rate 
for Event #4 (the "midnight flare" of 2002 April 15) has been determined independently by 
Torre et al.l (120121 ). who used a continuity equation analysis of the variation of the electron 



flux spectru m throughou t the s ource. The specific acceleration rate 7/(20 keV)= 11 x 10 3 s 1 



obtained by iTorre et al.l (120 121 ) is consistent with the value of 7.3 x 10 3 s 1 deduced here. 



The observationally-deduced value rj ~ 10 -2 s _1 implies that all available electrons 
would be energized and ejected towards the footpoints within a few hundred seconds. This 
result has significant implications for supply of electrons to the acceleration region, cur- 
rent closure, and the global electro dynamic environm ent in which electron acceleration and 
propagation occur (see, e.g.. lEmslie &: Henouxlll995l ). 



The values of the filli ng factor f deduced h e rein are broadly consistent with stochastic 



acceleration models (e.g., iPetrosian Sz Liull2004l ; iBian et al .1120 111 . Bian et al 2012, ApJ, in 
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press) which generally involve a near-homogeneous distribution of scattering centers. In ad- 
dition, the deduced values of the specific acceleration rates n are also broadly consistent with 
such models. For example, in their stu dy of electron (and proton) acceleration in a turbulent 
magnetohydrodynamic wave cascade, iMiller et al.l (119961 . their Figures 6, 7, 9, 10 and 12) 
derive values of the volumetric electron acceleration rate ~ (1.5 — 4) x 10 8 cm -3 s" 1 above 
20 ke V, with the exact va lue dependent on the assumptions of the various models considered. 
In the IMiller et al.l ( 119961 ) model, the background number density is n 



10 10 cm 3 , so that 



V 



1.5 - 4) x 10~ 2 s" 1 above 20 keV. 



On the other hand, values of the filling factor / close to unity pose significant chal- 
lenges for particle acceleration models that involve highly localized geometries, such as super- 



Dreicer acceleration in thin current sheets (see, e.g.. lLitvinenko fc Craigil2000tlTurkmani et al 



20061) . Turning to the specific acceleration rate in such an acceleration scenario. iHeerikhuisen et al 
( 120021 ) find (their equation [ 3.17]) a rate of proton ac celeration dj\f p /dt ~ 2 x 10 37 s" 1 , 
where £ (rj in the notation of IHeerikhuisen et al.ll2002l ) is the Lundquist nu mber, the ratio of 



the di ffusive to advective terms in the magnetic diffusion equation. Following IHeerikhuisen et al. 
(120021 ) . we take ( = 10 -8 , giving a proton acceleration rate dM p /dt ~ 2 x 10 33 s" 1 , and we take 
the number of protons available for acceleration as M p ~ en L 3 , where e ~ 0.4 y/C, — 4 x 10~ 5 
is the angle at the magnetic X-type neutral point at the origin of the acceleration region. 
With an ambient density n ~ 10 11 cm -3 and a longitudinal acceleration region extent 
L ~ 10 9 cm, M p ~ 4 x 10 33 and so rj ~ 2 s . Although this value is much greater than the 
values of n deduced here (it corresponds to the acc eleration of all amb i ent pa rticles in less 
than a second), it must be again stressed that the IHeerikhuisen et al.l (120021 ) model refers 
to highly efficient acceleration of a relatively small number of protons in a very localized 
geometry. We encourage calculations of specific acceleration rates for electrons in such a 
model, and indeed in all theoretical particle acceleration scenarios. 
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